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ABSTRACT 


A system for the analysis of geomagnetic and gradient 
noise is described. The nominal 176 kHz Larmor frequency 
output of an optically pumped cesium magnetometer is mixed 
to a frequency in the 1.2 kHz range and FM tape recorded. 
The recorded output is applied to a periodmeter, the digital 
output of which is input to the XDS-9300 computer for 
analysis and display. The absolute accuracy of the period 
measurement is ./5 usec with a resolution of .25 usec. 

The RMS noise of the tape recorder was determined to be 

.667 usec or in the operating frequency range less than 

1 Hz or .28 gamma. Externally introduced field disturbances 
as small as .5 gamma were recorded and displayed, and the 
response time of the magnetometer determined to be less 


than 2.5 msec. 
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I. INTRODUCTION 


With a continuing interest in magnetic anomaly detec- 
tion the Naval Postgraduate School has studied geomagnetic 
noise, first with the proton precession magnetometer and 
later with the optically pumped rubidium type. In 1972 
the school obtained an optically pumped cesium magnetometer 
with accessories to convert it to a gradiometer. In order 
to take full advantage of the inherent high frequency 
response of this magnetometer a data collection system 
was developed to observe geomagnetic and gradient noise. 
Figure 1.1 is a block diagram of the system. 

A brief introduction to the earth's magnetic field 
and other measuring instruments is presented first followed 


by a description of the system developed. 


I 





HEWLETT - 


3960 FM 
TAPE RECORDER 


MAGNETOMETER 


PERIODMETER 
AND 
INTERFACE 





XDS -9300 
COMPUTER 


AGT-10 
GRAPHICS 


Dalezbbers) FIR Ral System Block Diagram 





II. THE MEASUREMENT FIELD 


A, MAIN FIELD 

World magnetic surveys both at ground level and more 
recently by satellite indicate that the principal source 
of the main geomagnetic field is beneath the earth's crust. 
Approximately 90% of the main field exhibits characteristics 
roughly equivalent to the field of a terrestrially centered 


short bar magnet or dipole inclined at 11.5° to the axis of 


G 


rotation. See Figure 2.1. 






Figure 2.1 Dipole Appearance of Geomagnetic Field 
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The field strength is approximately 60 kilogammas in 
the polar regions, 30 kilogammas in the equatorial region, 
and 50 kilogammas in the neighborhood of 37° North latitude 
on the west coast of the United States. Figure 2.2 shows 
a world chart of the total intensity of the surface field.* 
The main field has a slow secular change of about .1% per 


year and it is a regional phenomenon which may reach .3% 


in certain regions. 


B. UNITS 


The customary unit of measurement is the gamma, regarded 
interchangeably as a unit of induction (gauss) or of magnetic 


intensity (oersted). 


1 gamma is 107? gauss (107? tesla) 
or 


1 gamma is 107? oersted (1072/47 ampere turn/meter) 


C. ELEMENTS 


The various elements of the magnetic field are shown 
in Figure 2.3. They are defined as: 


X North South Component 


lair Force Cambridge Research Laboratories, Handbook of 


Geophysics and Space Environments, p. 11-11, 1965. 
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World Magnetic Chart 





Y East West Component 

Z Vertical Intensity 

D Declination 

H Intensity of the Horizontal Component 
F Total Intensity 


I. Inclination or Dip 


ZENITH (-) 


“s NORTH (+) 





Figure 2.3 Elements of the Field 


In terms of the more familiar bar magnet: A magnet 
perfectly free to turn in any direction would show the 
vector direction of the field. It is almost impossible 


to suspend a magnet like that so two are used. One with 


M7 





a vertical axis (compass) and one with a horizontal axis 
(dip needle). The compass shows the direction of the 
horizontal component H, and when properly oriented the dip 


needle shows the inclination or dip. 


D. TIME VARIATIONS 

The time variations of the geomagnetic field can be 
broadly categorized into quiet variation fields, disturbed 
variation fields, and micropulsations. 

1. Quiet Variation Fields 

Geomagnetic records occasionally show a very smooth 

trace that clearly indicates a pattern of daily variation 
with respect to solar local time. This pattern is named 
"solar quiet day variation” and denoted Sq. There is also 
a lunar daily variation, denoted L, but it is much weaker 
and more variable. Both the Sq and L variations result 
from the dynamo action of the atmosphere, a complex phenom- 
enon of convective motions caused by the solar heating of 
the upper atmosphere under the influence of the earth's 
magnetic field. These convective motions induce ionospheric 
currents and electric currents in the earth itself, both 
of which cause magnetic variations. The effects of electric 


current induction within the earth are sometimes referred 
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to as geological disturbances because the amount of induc- 
tion is dependent on the composition of the earth's crust. 
For example, the intensity of induced currents in the ocean 
is higher than those in the continental crustal layers. 

Figure 2.4 shows the Sq variation recorded in 
Monterey during a quiet period. 

2. Disturbed Variation Fields 

Any fluctuation other than the quiet day solar and 
lunar variations can be called a magnetic disturbance. 
Other than disturbances directly attributable to a solar 
flare or fluctuations produced by upper atmospheric irreg- 
ular motions, all are caused originally by disturbed solar 
plasma and may therefore be generally classified as one 
family. The term magnetic storm is reserved to apply to 
a relatively severe, long lasting disturbance with recog- 
nizable features, but all other irregular variations 
deviating from the Sq and L variations are termed simply 
magnetic disturbances. They have durations lasting from 
100 hours to several minutes. Figure 2.5 shows a magnetic 
disturbance recorded in Monterey that immediately followed 


the Sq variations shown in Figure 2.4. 
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3. Geomagnetic Micropulsations 


Geomagnetic micropulsations are field variations 
with periods ranging from .2 seconds to approximately 10 
minutes. A terminology recommended by the International 
Association of Geomagnetism and Aeronomy is pce for pulsa- 
tion continuous and pi for pulsation irregular. The period 
range for these are: 
jee SL -27-) sec 
pe 2 5-10 sec 
peo 10-45 sec 
pe 4 45-150 sec 
pews 150-600 sec 
jpalee il 1-40 sec 
pi2 40-150 sec 
A general picture of the geomagnetic field in the 
lower frequencies is shown in Figure 2.6. 
Pe 1 pulsations have very small amplitude, between 
.Ol and .1 gamma and have a sinusoidal and beating form. 
The pulsations appear in bursts and usually last approxi- 
mately 30 minutes. They have been identified by their 
appearance and because there is a direct relationship 
between their reoccurrence period and mid-frequency. No 


seasonal or latitude dependence has been identified but 
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Figure 2.6 Geomagnetic Micropulsation Spectrum 


there is a diurnal variation in the mid-frequency events, 
high frequency events occur in post midnight hours, and 
low frequency events in the early afternoon. Figure 2./7 
is a record of a pe 1 micropulsation observed in Monterey. 
Figure 2.8 shows the typical form of the pe 2,3 
pulsations. They also exhibit a sinusoidal and beating 
character similar to the pe 1 type but their amplitude is 
slightly greater, usually .5 gamma. The duration of the 
signals seems to be related only to the sensitivity of the 
detector, which has lead to the hypothesis that when no 


Signal is observed the pulsation is not completely absent 
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1.8 gammas full scale; 45 second record 
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but rather too weak to be detected. Pc 2,3 events occur 
mainly during daylight hours and have a maximum near noon. 
The average period undergoes a diurnal variation, being 

a maximum at midday. They also have a 27 day (solar rota- 
tion) periodicity and a slight seasonal variation. 

As seen in the typical pe 4 record shown in Figure 
2.9, these micropulsations slowly rise to a maximum of 
between 5 and 20 gammas and then symmetrically decay. The 
average signal lasts about one hour. In mid-latitudes a 
daytime occurrence pattern can be observed but the period 
is shorter in the daytime than at night, the opposite of 
the pe 2,3. 

The longest period variations, pce 5, display ampli- 
tudes from 10 to hundreds of gammas. Their form is again 
Sinusoidal and they may last from several minutes to several 
hours. The largest amplitudes are observed in the auroral 
zone and the occurrences have a direct relationship with 
auroral zone phenomena. 

The class of micropulsations in the 2 to .005 Hz 
range with irregular form and a close relationship to upper 
atmospheric and auroral zone phenomena have been designated 
pi1l,2. Their amplitudes vary greatly, the larger sometimes 
being called micropulsation storms. Figure 2.10 shows the 


irregular appearance. 
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Figure 2.9 
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III. MEASURING INSTRUMENTS 


Today a wide range of devices are available for measur- 
ing the magnetic field, ranging from the extremely old 
magnetized needle to the newer superconducting quantum 
types. The major considerations in the selection of a 
device in a non-observatory application are sensitivity, 
accuracy, frequency response, mobility, and reliability. 
Three types of magnetometers that have developed with these 
considerations as guidelines are the fluxgate, optically 
pumped resonance, and superconducting quantum. They will 
be described and a comparison of the different varieties 


and other types is given in Table I. 


A. FLUXGATE 

A fluxgate magnetometer, sometimes referred to as a 
saturable core or saturable inductor, measures magnetic 
fields by utilizing the variation of permeability with 
intensity of magnetization in a ferromagnetic core. An 
easily saturable ferromagnetic core of high permeability 
is driven cyclically into saturation by an applied magnetic 
field. In the absence of an external field, which is 


usually de or very low ac and small with respect to the 
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applied field, the induced voltage in the sense windings 
will be symmetrical, i.e., contain only odd harmonics of 
the fundamental driving current that produced the applied 
field. In the presence of an external field, the induced 
voltage becomes asymmetrical, the asymmetry and the corre- 
sponding even harmonics being proportional to the amplitude 
of the component of the external field parallel to the axis 
of the core. 

The fluxgate magnetometer has been used widely in 
various applications where sensitivity, stability, and 
absolute accuracy has been sacrificed for the features of 
simplicity, low cost, low power, and ruggedness. Since it 
is a directional device, measuring only the component 
parallel to the axis of the core, a triaxial arrangement 


of three cores is necessary for total field measurements. 


B. OPTICALLY PUMPED RESONANCE 

The optically pumped resonance magnetometers (chiefly 
Rubidium, Cesium, and Helium) measure the total intensity 
of the magnetic field by an extremely accurate detection of 
the Zeeman effect. Due to the Zeeman effect the energy 
levels associated with a valence electron become split into 


sublevels whose separations are dependent on the total 





intensity of the magnetic field. The detection of this 
splitting is based on the principles of optical pumping 
and monitoring. Optical pumping is described in a tutorial 
manner in Reference 4 and for the Cesium atom in Appendix A. 
Figure 3.1 shows a simplified block diagram of a typical 
pumping apparatus. A light source producing an intense col- 
limated beam of light is first filtered and circularly 
polarized and then passed through a resonance cell containing 
the vapor to be pumped. As the light is absorbed and emitted 
by the sample atoms, optical pumping proceeds and results 
in a steady state non-equilibrium population distribution 
over the magnetic states of one of the ground levels. In 


this pumped condition the sample is prepared for resonance 


RESONANCE 
LENS FILTER CELL LENS 
PHOTO 
CELL 
LAMP ) 
CIRCULAR 
POLARIZER CULL 


Figure 3.1 Typical Pumping Apparatus 
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PmemlseChiaracter12ea Dy 42 Net magnetization M = M, along 
the ambient magnetic field direction. Resonance transitions 
between the Zeeman states of the pumped level are induced 
by subjecting the sample to an RF field of frequency f, 
sufficiently close to the resonant frequency £o and with 
component H, perpendicular to the ambient field Hy = Hz. 
The RF drive field introduces Zeeman coherence, and accord- 
ing to Bloch's equations the steady state behavior of the 
magnetization M changes and is no longer parallel to Ho but 
now has a transverse component M, precessing around Ho at 
the frequency of H,- Consequently the M, component along 
Hp decreases, and the light transmission also. This varia-~ 
tion of light intensity with the applied RF, being a maximum 
at resonance, can be used as a detection method, or the 
photocell can be made to observe the signal corresponding 
to the M, component, which is precessing or rotating at the 
same RF applied which induced the resonance. The photocell 
output is then at the resonant frequency. Appendix B 
details the two detection systems commonly used. The 
resonant frequency is proportional through a constant to 
the magnetic field (4.66 Hz/gamma for Rb; 3.499 Hz/gamma 


for Cs; 28 Hz/gamma for He) and thus the field measurement. 


Zo 





Figures 2.4 through 2.10 were all recorded using an 
optically pumped magnetometer and the appearance of the 
micropulsations differ somewhat from those one would find 
in a reference handbook. This illustrates an undesirable 
feature of this magnetometer when used in micropulsation 
work. It is equally sensitive to all frequencies and the 
larger amplitude, lower frequency changes such as the di- 
urnal variation continually drive the highly sensitive 
instrument offscale. The result is a sloping and discon- 


mmuous OUEpPUt that is difficult to interpret. 


C. CRYOGENIC SQUID 

A cryogenic SQUID (Superconducting Quantum Interference 
Device) magnetometer consists of three basic parts: 1) a 
Superconducting ring embodying a Josephson junction; 2) an 
externally generated HF field used to provide a small field 
variation superimposed on the ambient dc field and therefore 
to control the amount of flux through the ring; 3) an elec- 
tronic detection system to observe the impedance of the 
coil and make the output response a linear function of the 
ambient field. 

The operation utilizes the principle that if the flux 


through the ring is varied on either side of a minimum or 
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maximum energy value, associated with the Josephson coupling 
energy, the SQUID will generate a field which will respec- 
tively oppose or add to the variation. The energy contribu- 
tion of the SQUID aiding or opposing the applied field 
variation according to whether the ambient dc field amounts 
to a minimum or maximum energy ei constitutes a para- 
metric oreo of the inductance measured at the coil 
terminals. In other words, the SQUID alternately behaves 
diamagnetically and paramagnetically and will decrease and 
increase the inductive reactance. 

Through the use of a detector to monitor the coil imped- 
ance one can observe the variation of impedance with ambient 
magnetic field. Because the impedance variation is periodic 
a phase locked feedback loop is normally included to obtain 
a response linearly related to the field. 

Because SQUID compares the energy values of the ambient 
and weak superimposed fields, it measures only a change in 


field, not the total intensity. 
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D. GRADIOMETER 
A gradiometer, in theory, would measure the space 


gradients of the magnetic field, that is: 


Gee ot 4; 3x 
Ox ca OZ 
ay, ay , ay 
Ox oy OZ 
az, 22, 92 
ox Oy OZ 


Only five of these are independent because X, Y, and Z 

are the derivatives of a potential and because of Laplace's 
equation. In practice gradiometers are constructed by 
taking the finite difference of two magnetometers, there- 
fore cancelling the ambient field that is approximately 
equal at both sensors. A local change or anomaly that is 
closer to one of the two sensors is not cancelled and 
therefore recorded. 

The disadvantage of the gradiometer is that with respect 
to a single sensor magnetometer most magnetic anomalies fall 
off as the inverse cube, l/r? where r is the distance from 
the anomaly to the sensor, but the gradiometer anomaly falls 


off as the inverse fourth power, l/r where r is the 


a2 





distance from the anomaly to the center of the sensor 
configuration. 

In Reference 10, Langan, Lynch, and Slack reported a 
gradiometer used in geophysical prospecting that measured 
the vertical gradient A H/az where H was the magnitude of 
the anomaly and z the distance to the anomaly. It was 
constructed of two optically pumped sensors vertically 
separated and towed by a helicopter. 

Both the optically pumped and SQUID magnetometers are 
easily converted to gradiometers because of the form of 
their output, but they will measure different gradients. 
The optically pumped gradiometer, constructed of total 
field magnetometer sensors will measure the total field 
gradient AF/Az. A SQUID gradiometer, constructed of 
devices sensitive only to a field change will measure a 


field component gradient AZ/ Az, for example. 
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IV. SYSTEM DESCRIPTION 


The overall system can be broken into three subsystems: 
1) the data collection subsystem; 2) the periodmeter and 


~ 


interface subsystem; 3) the processing and display subsystem. 


A. DATA COLLECTION SUBSYSTEM 

The magnetometer used was the Varian V-4938 Cesium 
system. It used an optically pumped cesium sensor of the 
auto-oscillating type. A block diagram of the Varian sys- 


tem is shown in Figure 4.1. 


SENSOR 







LARMOR 


FREQUENCY 
~176 kHz 


200 FT 


READOUT UNIT 


Figure 4.1 Varian Magnetometer 


The sensor head produces the RF signal (Larmor frequency) 
directly proportional to the total magnetic field intensity. 
This signal together with the input power for the sensor is 


conducted along a 200 foot coaxial cable to the main readout 


5) 9. 


unit. The readout unit converts the Larmor frequency to a 
highly stable de voltage which drives the pen recorder. 
The system was located in a wooded area approximately one 
mile from the Naval Postgraduate School, which provided a 
near noise free environment. 

A more complete block diagram of the V-4938 system is 
shown in Figure 4.2. The RF Larmor frequency, which was 
approximately 176 kHz, from the sensor is mixed in the read- 
out unit to the 100 Hz to 3 kHz range and then applied to 
a highly stable frequency to voltage converter. This dc 
voltage, which now represents the field value, is passed 
through biasing circuitry to interface it to the pen recorder. 

The system configuration used for data collection differed 
from the Varian recommended system of Figure 4.2. This was 
necessary because of the requirement to record the output 
on magnetic tape. The actual system used is shown in Figure 
4.3. This configuration was used because the mixer output 
of the readout unit was a square wave and its frequency 
content exceeded the 5 KHz range of the tape recorder. The 
coupler unit shown connects between the readout unit and 
the sensor to provide dc power to the sensor and it also 
contains a mixer, but its output is a sinusoid rather than 


a square wave and recorded nicely. The coupler units primary 
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Figure 4.2 Magnetometer Signal Flow 
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Figure 4.3 Magnetometer Data Collection Configuration 
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purpose is to convert the system to a gradiometer which is 


shown in Figure 4.4. 


7 SENSOR #1 ] SENSOR #2 
Bs 
FREQ. 
REF 


| te he COUPLER #2 
a | ite 


MIXER OUTPUT 


MIXER OUTPUT MAGNETOMETER 
GRADIOMETER 





Figure 4.4 Gradiometer Data Collection Configuration 


Again this is not the configuration recommended but 
this provided both a magnetometer and gradiometer output 
suitable for recording and still allowed the pen recorder 
to monitor the gradiometer signal. The bias applied to 


the one sensor was required to generate a difference 
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frequency between the sensors sufficient to drive the mixer 
circuitry. The sensors were located approximately 180 ft 
apart resulting in a Larmor frequency difference of only 
approximately 30 Hz, and the mixers required a difference 
of at least 100 Hz and optimally between 1 and 3 kHz. 
Figure 4.5 shows the readout units mixer output with an 
insufficient difference frequency and Figure 4.6 shows the 
coupler units sinusoidal mixer output. 

One other system characteristic imposed a limitation, 
that was the orientation sensitivity of the single cell 
sensor. For optimum performance the sensor head needed to 
be oriented both toward magnetic North and 45° to the total 
field vector. Figure 4.7 shows the non-metallic housing 
and Figure 4.8 the geometry. 

Figure 4.9 shows the system components and equipment 


used to monitor and record the data. 
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Figure 4.5 Distorted Readout Unit Mixer Output 





Figure 4.6 Coupler Unit Mixer Output 
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Figure 4.8 


Sensor Housing Geometry 
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Figure 4.9 


Pen Recorder, Coupler Units, and Tape Recorder 





B. PERIODMETER AND INTERFACE SUBSYSTEM 

The frequency response of the magnetometer when used 
With the readout unit is limited to 20 Hz by the discrimin- 
ator and even more by the pen recorder response. The 
theoretical frequency response of the sensor is limited 
only by the bandwidths of the amplifier and feedback loop, 
which are extremely wide. To take advantage of this the 
mixer output was recorded on an HP-3960 FM tape recorder. 
This signal still contained the high frequency components 
present in the Larmor signal. A periodmeter and interface 
unit were constructed to measure the period of successive 
cycles and interface this digital count into the XDS-9300 
computer, for further display on the AGT-10 graphics console. 

Figure 4.10 is a block diagram of the periodmeter. The 
data signal is passed through a simple RC high pass filter 
to reject any 60 Hz modulation that might be present and 
then input to a zero crossing detector. The falling edge 
of this is used to trigger a one-shot that inhibits the 
otherwise free running 4 MHz count. The one-shot pulse 
is inverted and delayed and the rising edge used to trigger 
flip flop type memories. Once the count is latched another 


one shot is used to clear the count and provide a signal 
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indicating the cycles completion. The cycle then repeats. 
The actual circuits and timing diagrams are presented in 
Appendix C. 

A direct access channel into the XDS-9300 was used to 
input the binary count directly into sequential memory loca- 
tions. Figure 4.11 is a block diagram of the periodmeter 
and interface. To implement the input the following sequence 


of instructions were used: 


EOM 032400 
SKS 032400 
BRU $-1 
el TU bai 


With the execution of the EOM instruction the contents of 
the control register are available at the output connector. 
This is decoded and toggles a flip flop to alert the period- 
meter that the XDS is ready to accept an input. This EOM 
alert and the periodmeter's completion signal are logically 
combined and send a ready signal to the XDS-9300. Upon 
receipt of this ready signal and the execution of the PIN 

the data is transferred and stored in the location specified 
in the PIN instruction. The SKS instruction is used to sense 
if both the computer and periodmeter are ready. If both the 


EOM alert and completion signal are not present, the next 
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instruction, branch unconditionally back to the SKS, is 
executed. This looping continues until both systems are 
ready. If both are ready the branch instruction is skipped. 
The ready signal is also used to gate the data through in- 
verters because it is inverted prior to being transferred 
into the XDS-9300. The actual circuit diagrams are pre- 


sented in Appendix C. 


fee DESPLAY SUBSYSTEM 
A block diagram of the data flow in the interface and 


display subsystems is shown in Figure 4.12. 


PERIOD - 





Figure 4.12 Data Flow From Periodmeter Through Display 


The physical block length of data input to the XDS-9300 
was 12,288 points. A mean and standard deviation of these 
points was found and values not within five standard 


deviations were discarded. This eliminated erroneous data 
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points caused by noise spikes: either in the computer, 
periodmeter, or on the record:ng. The block was then 
buffered out in groups of 4096 points and written on 
magnetic tape. To be displayed the data was buffered 
back in, again in groups of 4096 points, normalized, and 
Bheelaved on the AGT-10 graphics console. Both programs 
are self-explanatory and are presented in Appendix D. 
Signal flow ee for the two are shown in Figure 4.13 


and 4.14. 
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Figure 4.13 Signal Flow Diagram of Input Program 
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Figure 4.14 Signal Flow Diagram of Display Program 
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V. RESULTS AND CONCLUSIONS 


The system was tested in three phases: 1) laboratory 
signals; 2) tape recorded laboratory signals; 3) magnetom- 
eter signals. 

First, constant frequency signals from a laboratory 
generator were applied. The frequencies sampled and the 


resulting period measurement are listed in Table 2. 


FREQUENCY Hz MEAN PERIOD usec THEORETICAL PERIOD psec 


1000.0 999.25 1000.00 

1000.1 999.00 999.00 

1000.5 998.75 999.50 

1001.0 998.25 . 999.00 

1003.0 996.25 997.01 

1005.0 994.25 995.02 
Table 2 


The .75 usec difference between theoretical and actual is 
due to the clear and latch time of the periodmeter. These 
measurements were within +1 count or +25 wsec with a 4MHz 
time base. This established the resolution of the device 


as .l1 Hz or better. The outputs displayed for the sampled 
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frequencies of 1000.0 Hz, 1000.1 Hz, and 1001.0 Hz are 
shown in Figures 5.1 - 5.3 and a distinguishable difference 
is noticeable. To insure the system would follow a rapidly 
changing frequency a swept frequency signal was applied and 
the display shown in Figure 5.4. 

The second phase of testing was designed to determine 
the magnitude of the tape recorder flutter and noise. It 
involved recording the same frequencies used in phase 1 and 
comparing the results. The frequencies, mean period measure- 


ment, and standard deviation for each are listed in Table 3. 


FREQUENCY Hz MEAN PERIOD usec STANDARD DEVIATION uS ec 


1000.0 DoS je 5 ee 
KOOO =a 99,9500 . /00 
1000.5 D164 1 5 ~850 
1001.0 I9 5225 D0 
1003.0 BEN 5 72 - 650 
1005.0 9G 325 nodD 


The mean period, calculated for 12,288 periods, remained 
the same. The average standard deviation was .667 usec. 
This established the RMS noise level of the system when 


operated in the 1200 Hz range to be approximately 1 Hz or 
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Figure 5.2 Output Display of Direct Input 
1000.1 Hz Laboratory Signal 
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.27 gamma. The AGT-10 displays for the sampled frequencies 
of 1000.0 Hz and 1000.1 Hz are shown in Figures 5.5 and 5.6. 
The previously discernible difference is now buried in tape 
recorder noise. Figure 5./ shows the frequency transition 
from 1001.0 to 1003.0 Hz and this 2 Hz difference is easily 
distinguishable. 

The final phase of testing involved introducing a known 
magnetic disturbance and analyzing the magnetometer output. 
A 72 turn Helmholtz coil, in the configuration of Figure 5.8, 
was used to introduce the disturbance. The sensor was 
placed on the axis of the coil and the coil aligned with 
the total field vector. This insured the induced field of 
the coil would directly aid or oppose the ambient field 
and allowed changes of predetermined amounts. With the 
simplified geometry 
SRE tng 

Powe Tae 
where N = number of turns in the coil 
I = current in the coil 
R = radius of the coil 


Ee=aistance £rom coil tO sensor 


a was required for conversion from 
Os ampere-turns per meter to gammas. 


Sa 
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Figure 5.6 Output Display of Recorded 
1000.1 Hz Laboratory Signal 
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Figure 5.7 Output Display of Recorded1001.0 


and 1003.0 Hz Laboratory Signal 
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The coil utilized had a mean radius of 33 cm and was placed 
3.5 m from the sensor. Table 4 tabulates the various 
currents and observed versus predicted changes. Figures 
5.9 through 5.30 show the magnetometer pen recorder output 


and a selected time frame of the corresponding AGT-10 display. 


THEORETICAL MAGNET. THEORETICAL DISPLAY 

I ma OH AH ian fas 1 
400 45.964 46.0 160.874 160.893 
300 34.473 Doe 0 ZO Oo) 120.666 
240 UR SIMs. Jes: 96; 524: ee yeiss) ik 
200 DO Zoey 80.437 80.442 
160 18.386 loa, 64.350 64.354 
120 SiG eee, lige Pare: 43.262 43.264 
100 11.491 ie el 40.218 40.219 
90 10.342 54S Ste US 36.198 
80 oon al Shs ae Sy 6 lyfe oi 6 
70 8.044 8.4 2315159 28.154 
60 On.095 6.8 24.131 ZO LOZ 
4.0 4.596 4.2 Ie) (Orery 16.086 
20 22296 Jaa 8.044 8.045 
10 1.149 1D, il 4.022 4.023 
5 574 05 Zeon L ZOnLZ. 

Table 4 
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Figure 5.11 Pen Recorder Output for AH= 35 gamma 
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Figure 5.12 AGT-10 Display for MAH=35 gamma 


68 





a 
tt 


oe == == 525 =2=========== 













SSS SPs S=Es 

, aE SSeS 
: =: S=SS==! 

: Briss == ==! 
= 

= 

aa 

QO. 


TIME, SEC 


Figure 5.13 Pen Recorder Output for AH= 28 gamma 
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Figure 5.14 AGT-10 Display for 4H= 28 gamma 
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Figure 5.16 AGT-10 Display for AH= 18 gamma 
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Figure 5.17 Pen Recorder Output for AH= 11 gamma 
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Figure 5.18 AGT-10 Display for AH= 11 gamma 
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Figure 5.19 Pen Recorder Output for OQH= 8 gamma 
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Figure 5.20 AGT-10 Display for AH= 8 gamma 
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Figure 5.21 Pen Recorder Output for AH= 4 gamma 
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Figure 5.22 AGT-10 Display for AH= 4 gamma 
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Figure 5.23 Pen Recorder Output for AH= 2 gamma 
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Figure 5.24 AGT-10 Display for AH= 2 gamma 
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Figure 5.25 Pen Recorder Output for AH= 1 gamma 
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Figure 5.27 Pen Recorder Output for AH= .6 gamma 
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To obtain a precise measurement of the magnetometer's 
response time the AGT-10 display was expanded. Figure 5.29 
shows the rising edge and Figure 5.30 the falling edge. As 
seen the rising edge required approximately three periods 
or 2.5 msec to respond and the falling edge only one period 


or .8 msec. 
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Figure 5.29 Expanded View of Rising Edge 
of Magnetometer Response 
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Figure 5.30 Expanded View of Falling Edge 
of Magnetometer Response 
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The periodmeter and interface unit measured the period 
of the signal with an absolute accuracy of ./5 sec and 
properly input the value to the computer. The resolution 
of the device was .25 fA Sec using a 4 MHz time base and 
this could be increased by using a faster time base. The 
system was successfully operated with a 10 MHz clock fre- 
quency and resolution of .1l sec. The use of an increased 
clock frequency or averaging of periods to increase the 
periodmeters resolution would not improve the system sensi- 
tivity. The overall system sensitivity is limited by the 
flutter and noise of the tape recorder. Two methods were 
attempted to overcome this: 1) the use of a synchronized 
ayscillator constructed of two multivibrators synchronized 
with the zero crossings of the input signal; 2) recording 
a 1 kHz reference signal on a second channel of the 
recorder and multiplying it up for use as a time base. 

The pulse width instability of the multivibrators was 

greater than the flutter noise so the synchronized oscillator 
was unsatisfactory. A phase lock loop type frequency multi- 
plier was constructed to multiply the 1 kHz to 1 MHz but 

the delay time through the loop was too great to cancel the 
flutter for a period to period measurement. 

To achieve greater system sensitivity and resolution some 


means of flutter cancellation is necessary. 
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APPENDIX A 


OPTICAL PUMPING OF THE CESIUM ATOM 


Figure A-1l shows the energy levels of the cesium atom. 
Absorption of photons of light by the cesium valence elec- 
trons cause their energy to be increased from the “a ground 
state to the “Py state. Due to the light's polarization 
the magnitude of the quantum number m is also increased by 
one. An excited atom quickly emits its energy and returns 
to any of the ground state sublevels. The absorption and 
re-emission is repeated resulting in a shift to the m=+4 
sublevel. Since the “Py state contains no m level above 
+4, transitions out of the “Si, =+4 level are not possible. 
This is precisely optical pumping and under these conditions 


there can be no further absorption of light. 
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Figure A-1l Energy Levels of Cesium Atom 
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APPENDIX B 


OPTICAL DETECTION SYSTEMS 


The first method used to obtain a continuous record of 
the total field intensity used a control loop arrangement 
and was therefore called a Controlled Magnetometer. A 


block diagram of it is shown in Figure B-l. 
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Figure B-l Controlled Magnetometer 


A small frequency modulation is applied to the local RF 
oscillator to sweep it back and forth through resonance. 
Owing to this sweeping the signal collected at the photo- 


cell is also modulated at the same frequency, but its phase 
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with respect to the initial modulating signal depends on 
the position of the local RF with respect to the resonant 
frequency. The photocell output and modulating signal are 
phase detected and the output is used as an error signal 
in a control loop to maintain the local RF oscillator coinci- 
cent with the resonant frequency. This system does not place 
severe restrictions on the performance of the photocell or 
its amplifier, has high accuracy, good stability, and good 
low field capability. 

One attribute the controlled system does not have is a 
good response to transients. Figure B-2 shows the block 
diagram of a system that has the excellent response to 


transients necessary to detect and observe micropulsations. 
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Figure B-2 Self-Oscillating Magnetometer 
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In theory a second light beam is placed at right angle with 
respect to the first and the variation of the M, component 
of the magnetization is observed. Since this transverse 
component oscillates at the frequency of the signal injected 
and the amplitude is a maximum at resonance, the photocell 
is simply amplified, phase shifted, and re-injected to the 
coils. The system will then oscillate at the resonant 
frequency. In practice it is possible to use a single light 
beam inclined at 45° to the direction of the field to be 
measured. This places an orientation restriction on the 
sensor, but systems employing two and three resonance cells 
have been developed that overcome this. This system places 
severe restrictions on the amplifier and especially on the 
accuracy of the phase shift, which is required to maintain 
oscillation. It has a low absolute accuracy and poor long 
term stability but as previously mentioned, excellent 


response to transients. 
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APPENDIX C 


PERIODMETER AND INTERFACE CIRCUITS 





P 


TL IL 


7407 
74175 
74163 f 


Latch -E 


Inhibit 
Clock 
Clear 


Figure C-1 Counting Circuitry 
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Figure C-2 Input Circuitry Timing Diagram 
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